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We find general solutions to the generating-function equation c. 



F(z) 



where X is a 



complex number and F(z) is a convergent power series with F(0) 7^ 0. We then use these results to 
derive finite expressions containing only integers or simple fractions for partition functions and for 
Euler, Bernoulli, and Stirling numbers. 
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I. INTRODUCTION AND MAIN RESULT 

Generating functions are often used as a compact way to define special number sequences and functions as the 
coefficients in a power-series expansion of more elementary functions. One example of this is the partition function 



p(n), (sequence A000041 in OEIS [l[ 
function 



the number of partitions of n into positive integers. It has the generating 



9 =o 



CO 1 



1 + z + 2z 2 + 3z 3 + 5z 4 



k=l 



Some other examples of generating functions are those for the Bernoulli and the Euler numbers: 
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y^-zj = 



-1' 
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2q_ 2q _ 



cosh; 



(1) 



(2a) 



(2b) 



The solutions of generating-function equations can often be found from recursion relations. Some of them can also be 
calculated directly by using various expressions, such as Laplace's detcrminental formula for the Bernoulli numbers 
Q, which we write in the form: 



B n = nl 



1 

2! 



(«-!)! 



(n+l)l 



2! 



(3) 



Vella 4] has derived expressions for B n and for E n as sums over the partitions (or, alternatively, over the compositions) 
of n using a method based on the Faa di Bruno formula for the higher derivatives of composite functions. Concerning 
the partition function, Rademacher 0, Q derived an exact formula, an improvement over the Hardy-Ramanujan 
asymptotic formula, which however involves an infinite sum of rather complicated, non- integer terms. More recently, 
Brainier and Ono [6| have derived an explicit formula for the partition function as a finite sum of algebraic numbers 
that requires finding a sufficiently precise approximation to an auxiliary function. 

In the following theorem, we present methods for solving generating-function equations for the case where the 
right-hand side is expressible in terms of a convergent, non-zero-near-the-origin power series. We will then apply this 
theorem to solve for the partition function and other objects which have the appropriate generating functions. In the 
following, the notations 

<5n,<7iH \-q K 



<Zi > • • • > <1k J qi ! • • • gif! 
and 

fX\_ X(X — 1) ■ ■ ■ (X — N + 1) 



N 
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will denote, respectively, multinomial coefficients and generalized binomial coefficients. 
Our main result is: 

Theorem 1 Let F(z) be a holomorphic function in a neighborhood of the origin with the series expansion Y^Lo a q z9 > 
with a 7^ 0, and let X be a complex number. Then the coefficients c p X ^ in the generating-function equation 



x 



p=0 

are given by the two equivalent expressions: 



(I) 



E 

0<k\ ,...,k p <p 



x 



K 



^p,J2 mk m a Q 



X-K fci 



where K = k\ H + k p , and 



(II) 
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The proof of (II) requires some additional machinery and will be postponed to Section II. 

Proof of (I): Since (ao + a\z + ■ ■ ■ ) x — (a + a\z + • • • + a p z p ) x ~ z p+1 , the coefficient of the z v term is the same 

( X) 

in the two sums, so in determining the c p coefficient we need deal only with the finite sum. Using the generalized 
binomial theorem, we write 



x 



e«« 

\q=Q 



E 

oo 

E 

N=0 



X 

N J a 'o 



N 



X-N 



X 

N ) a 'o 



X-N 



E a ^ 9 

\q=l ) 

E 

0<fei,...,fep<7V 



N 

k\ , . . . , k p 



■a p p z ^ mk " 



In the second line we've applied a multinomial expansion to (X)g=i a q z q ) N . We now write z^ mk " 
and interchange the order of the sums: 



x 



E°« 

\q=0 



E* s E E 

s=0 0<ti,...,fc p <sJV=0 



N 



fa, 



, k p 



c X-N k^ 



(4) 
(5) 



After the interchange, the fc g 's are fixed in the sum over N and only the TV = k\ + ■ ■ ■ + k p (= K) term contributes. 
We are interested in the coefficient of z p in this sum, which is thus 



C (X) = 



E 

0<k± ,...,k p <p 



K 

k\, . . . , k p 



X X-K fei 

o p ,J2mk m a a 1 



■a k /. 



(6) 



As stated above, this is also the coefficient of z p in the infinite sum. QED. 

The Kroneker delta in eq.(6) restricts the sums over the fc's to a sum over the partitions of p. The multinomial 
coefficient in this formula counts the number of unique ways the parts can be ordered. As a sum instead over 
compositions, eq.(6) becomes 



"P,<?iH ^<1k a a <l 



<i x> = ±{ x K ) e • 

K=0 V 7 l<qi,...,q K <P 

As a corollary, we have the multinomial identity: 
Corollary 1 For fixed p, 

E(-d* E L p , qK )=(-v p 



(7) 
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Proof: We set X = — 1 and a q — l/ql, so that (J2 a q zq ) X = e 2 and so a 



(-i) 



-l) p /p\. Then, from eq, (7), 



p! 



EM)" E 



l<gi,...,gjf <p 



u p,gH hgK 



(8) 



Multiplying both sides by p\ gives the result. 
QED 

We now apply Part (I) of Theorem 1 to eq.(l) and use Euler's pentagonal theorem: 



JJ(1 - z k ) = 1 - z - z 2 + z 5 + z 7 - z 12 - ■ ■ ■ . 



(9) 



k=l 

where the exponents 0, lj2, 5, 7, 12, . . . are generalized pentagonal numbers (sequence A001318): q m = (3m 2 
m)/2, m = 0, ±1, ±2, • • • [2j. p(n) is then equal to a sum over the pentagonal partitions of n: 



p(n) 



E (-D 



fcl, K2, K5, . . . , k, 



>n,J2 k qm q„ 



(10) 



where qu is the largest GPN < n. Eq.(10) thus expresses pin) as a finite sum of integers. The number of terms in 
the sum is the number of partitions of n into generalized pentagonal numbers, (sequence A095699). For example, 9 
has 10 pentagonal partitions, (9 = 7 + 2 = 7 + l + l = -- - = l + -- - + l), and p(9) is 



p(9) = - 



2 

0,1,0,1 
5 

1,4,0,0 



3 

2,0,0,1 
6 

3,3,0,0 



3 

0,2,1,0 
7 

5,2,0,0 
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2,1,1,0 
8 

7,1,0,0 



4,0,1,0 
9 

9,0,0,0 I • 



From the as-yet-to-be-proven Part II of the theorem, p(n) is also expressible as a (k + 1) x (k + 1) determinant, for 
any integer k > n; 



p(n) 



10 

-110 

-1-1 1 

0-1-1 1 



dk-i dk-2 dk-3 dk-4 
dk dk-i dk-2 dk-3 
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1 : 
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(11) 



where, for < q < k, 
d q = (sequence A010815) 



(-l) m if q = q m , rn = 0,±l,±2, 
otherwise, 



This reduces to the form stated in the theorem by successive expansions by minors along the top row. Form (11) will 
prove more useful in some of the following discussions, but it can be reduced (by expansions in minors) to the n x n 
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determinant, 
GPN's 
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II. MATRIX FORMALISM 



Matrices which are constant along all diagonals are Toeplitz matrices. We will be concerned in this section with 
lower-triangular Toeplitz (LTT) matrices. A nondegenerate LTT matrix is one with non-zero diagonal elements; if 
the diagonal elements are equal to 1, it is then a unit LTT matrix. 

Infinite-dimensional LTT matrices have the form 



A = 



a 

a 2 ai a 







v 



(13) 



The determinant of a fc-dimensional lower-triangular matrix is the product of it's diagonal elements; if the matrix is 
also Toeplitz, then its determinant is Oq. 

We define the infinite-dimensional lower shift matrix J as 



/ 







v 



\ 



The pth power of J has elements (J p )i 



(14) 



,•; these matrices obey the relations J p J q = J q J p = J p+q . Any 



infinite-dimensional LTT matrix can be expanded out in non-negative powers of J: 

U \ 





Oi a 
a 2 a\ a 



V 



q=0 



(15) 
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The product of two infinite-dimensional LTT matrices is 



AB = J2 a r jq J2 bsJS = J2 JP J2 a i b P-i = BA 

q=0 s=0 p=0 q=0 



(16) 



A similar expression can be written for finite-dimensional LTT matrices, with the J's replaced by finite lower shift 
matrices. LTT matrices, whether finite- or infinite-dimensional, thus commute with one another. 

Lemma 1 Let A be an infinite- dimensional nondegenerate LTT matrix: 

oo 

A = ^2a q J q , a ^0. 

q=0 

Then the inverse of A is an LTT matrix with coefficients 



b p = 



ai 





a 











«o 
ai 



for any integer k > p. 
Proof: We have 

oo oo oo p 

q=0 s=0 p=0 q=0 

But 



(17) 



q=0 



,fc+l 



a 
ai 





a 



a 



a 

Oi 



1 if p = 0, 
if p > 0, 



(18) 



and therefore J]^ 6 S J s = A" 1 . QED 

It is clear that this proof still holds if the J's are replaced by finite-dimensional lower shift matrices and the sums 
over q and s are finite. Thus, the coefficients of the inverse of a finite fc-dimensional (nondegenerate) LTT matrix are 
given by the same formula, and we have the result: 

Lemma 2 The inverse of a k-dimensional LTT matrix is equal to the k-dimensional truncation of the inverse of the 
corresponding infinite- dimensional LTT matrix. 

I.e., 



/ bo 

bi b 

b 2 h b 

V : 

iff 

( bo 

h b 

\ bk-i bk-2 ■ ■ ■ 
for all k. 











a 

ai a 
a 2 ai a 

V : 



a\ ao 



bo J 







N- 1 



y afe_i afe_2 • • • 
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Lemma 3 If A = J2^Lo a pJ p * s an infinite- dimensional nondegenerate LTT matrix and X is a complex number, then 
A raised to the power X is 



^ = E^ S E 

s=0 0<ki,...,k s <s 



X 
K 



K 

k\, . . . , k s 



A ^ n X ~ K n kl n k > 

O s ,J2mk m a a l a s 



Proof: As before, we first consider the finite-sum case. Let z be some nonzero complex number. We then define Z 
as the infinite row vector 



Z = (l, z, 



•)• 



(19) 



Multiplying Z by J on the right, we have ZJ = zZ. Therefore, Z J" = z^Z and, by linearity, ZQT c q Ji) = c g z«)Z. 
Then, if F(z) is a holomorphic, nonzero function of z in a neighborhood U of the origin, ZF(J) = F(z)Z for all 
z € U, where F(J) is the LTT matrix obtained by replacing z by J in the Taylor series expansion of F(z). Since 
(Y^iq—o a p zP ) X > a o 7^ is such a function, we can write, 



x 



x 



z Emm = Em 



\<?=0 



\q=0 



= E^ S E 

s__0 0<fei,...,fc p <s 
oo 

^E^ s E 

s=0 0<fei,...,fe p <s 



X 
K 



K 

fcl j • • • j fcp 



<5«,£n»k m a a i • • • V z 



A' / V k\ , . . . , kp 



A „ n X ~ K n kl n k P 

0s,J2mk m a Q a l ' ' ' a p 



z is otherwise arbitrary, so we have 

(p \ X oo 

e^m =e^ s e 
<j=0 / s=0 0<fei,...,fe p <s 



X 
K 



K 

k\ , . . . , kp 



A „ n X ~ K n k ^ r, k P 

s ,J2mk m a l " " " fl p 



For fixed s < p, k s+ i = ■ ■ 
( p \ x p-i 

e^m -e^ e 



\q=0 



s=0 0<fei,...,fe s <s 
oo 

+ E^ S E 

s=p 0<fei,...,fe p <s 



X A' 
K \ Id 



, . . . , k s 



if / V k!,...,k p 



c X-K fei fc s 



<Js,*jTmk m a a l ' ' ' a p ■ 



The equality in the lemma is then demonstrated by letting p — > oo in this equation. QED 
We now have the necessary machinery to prove the second part of Theorem 1. 



(20) 



(21) 



k p = in the sums over the fc's as a result of the restriction s = ki + ■ ■ ■ + pk p : 



(22) 



Proof of Thm 1,(11): From the proof of Lemma 3, the coefficients in the expansion of (^2a p J p ) X in powers 



of J are the same as in the expansion of (J2 a p z p ) X in powers of z; i.e. 



E c ^ v - E m iff E c o jq = E a * jq 



9=0 



\q=0 



\q=0 / 9=0 
We can therefore use Lemma 1, written in the form 

/ ao \ / 

ai a 



b p - v+l 



a p -i 

V a p a p-i ■ ■ ■ 



a 

Oi a J 







(23) 



1 \ 





(24) 



7 



(taking k — p), with the replacement A — a p J p —tA= (J2 a p J p )~ X (and thus ao —> l/a x ), to solve for c p X ^ as a 
determinant in the form (II) in the theorem. QED. 
In a straightforward fashion, eq. (23) generalizes to 

oo / oo \ Xl / oo \ X " 

E^ P = £<H ••• E4"^ 



p=0 



Vp=0 



\P=0 



iff £ Cp J? = ^ fl Wj" 

p=0 Vp=0 / 

The c p coefficients are then: 



\p=0 y 



(25) 



P+i 



/a (1) 



(1) Jl) 



ai a, 



o 







+ 



( 





(«) («) 







(«) _(n) 
P "p-l 



4 n) / 



V 



-(a«)- Xl 



(26) 



Gradshteyn and Ryzhik [8j give an equivalent expression for these coefficients for the case n = 2, X\ = 1, X2 = — 1. 
Returning to the partition function: MacMahon's recurrence relation Q, 



p(n) — p(rt — 1) — p(n — 2) + p(n — 5) + p(n — 7) — p(n — 12) — p(n — 15) + • • • = 0, 



(27) 



follows directly from expression (11); setting k = n in that equation, the sum on the right in (27) is equal to the 
determinant 



1 

-1 1 

-1 -1 

-1 







1 

-1 1 







1 d„_ 



d n d n -\ d n _2 d n . 



-1 



which is zero since the first and last columns are equal. 

Expressed in terms of the J matrices, expression (12) for p(n) is 



p(n) 



dot 



q m <n+l 



_|_ (— l) m +! ^ j(m-l)(3m+2)/2 _|_ j(m+l)(3m-2)/2 ^ 



m>0 



(28) 



where J T is the transpose of J and the notation [ } nxn means the n x n truncation of an infinite-dimensional matrix. 
And, by relation (23), we have the compact matrix equivalent of p(n)'s generating-function equation: 

••• o\ _1 

110 
1-11 



( p(o) 

p(l) P (0) 
p(2) p(l) p(0) 



\ 





/ 1 



(29) 



V p(k) P (k - 1) p(k-2) ■■■ p(0) ) \d k 4-i 4-2 • • • 1 / 
The generating function for the number of partitions in which no part occurs more than D times is Q 



n 

k=l 



1 _ X {D+X)k 



1 - X 



D+l 



r 2(D+l) 



X 



5(D+1) 



X 



7{D+1) 



1 — x — X 2 + X 5 + X 7 



(30) 
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We have then from (26), with X\ = = 1, 



Pd(ti) 



where 



to = 



1 

-1 1 

d n -l d n ^2 
d>n d n —\ 



t 

h 

1 tn-l 
-1 t« 



(31) 



(-l) m if q = (D + l)q m , m = 0, ±1, ±2, 



(32) 



otherwise. 

By expanding this determinant by minors along the last column, and using the expression (11), we get the relation 
PD (n) = p(n) + ]T(-l)>(n - (D + l)q n ), (33) 

m 

which is a generalization of (27), the D = case. If we now take D = 1, pAn) = </(n) = the number of partitions of 
n into distinct integers (A000009). This is also, from a result due to Euler the number of partitions of n into odd 
integers, and so we have for the odd partition function 



q(n) = p{n) - p(n - 2) - p(n — 4) + p(n — 10) H h (-l) m p(n - 2p m ) 



(34) 



III. APPLICATION TO BERNOULLI, EULER AND STIRLING NUMBERS 

Vella's expression for the nth Bernoulli number, (eq.(a) in his Theorem 11) is, in my notation, 



Bn = y { -^- y 

^ l + K ^ 

K=0 l< qi ,...,q K <n 



n 



Applying Part I of Theorem 1 to the generating function (2a) for Bernoulli numbers, we have 

k n 



K 

k\ , . . . , fc„ 



^ \ ki / * \ k2 



0<ki,— ,k n <n 

As a sum over compositions, this expression becomes 



2! 



3! 



-1 



(n + 1)! 



Bn = EM)" £ 



K=0 



l<9i,. ..,q K <n 



(qi + 1) • • • (qk + 1) \qx,...,q K 



(35) 



(36a) 



(36b) 



to be compared to Vella's result above. 

Another expression for the (even-numbered) Bernoulli numbers is obtained from the generating function 

^ (2q)\ sinhz' 



(37) 



From this we have 
(2p)l 



B 2p — 



£ 



-1 



(2p+l)! 



^EH) K E 



2p 



(38a) 
(38b) 
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while from the generating function for the Euler numbers (A000364), eq.(2b), we obtain the two expressions, 

«*-<*»» £ ( J^-^)* 1 (=wT-(m^ • <39a) 



0<fci,...,fc p <p 

E(-!) K E f 2 2 \ )■ ( 39b ) 



i?2p and Ei v can thus both be expressed as sums over the even compositions of 2p. Since there are no even 
compositions of odd numbers, these expressions can be extended to include all of the odd-numbered numbers except 
for B\. We then have, more generally, 

B - - ^|-^j„s»(^:J' (40) 

(n + 1), 

Ln/2j 



f £ 1 

if=0 l<«i,...,gjf<Ln/2J 

This expression for _E„, in a different notation, was previously derived by Vella, (eq. (c) in Thm. 11 4]). E n can 
however also be expressed as a sum over the odd partitions/compositions of n — 1: 

Proposition 1 For n >1, 



L»/2J 

En = E E 

N = l l<gi,.... 92 N-l<|n/2j 



n- 1 

2qi - 1, . . . , 2q 2 AT-i - 1 



Proof. The sum is over all compositions of n — 1 that contain an odd number of odd parts, which is an empty set 
if n — 1 is even. Therefore E n = for odd n, and we only have to prove this equation for even n. The equality is 
trivially true for n = 2. We will show that both Ei v and 

iV=l l^,..,^-!^, V ^ - 1> ■ • ' ' 2 «2JV-1 - 1 / 

satisfy the same recursion relation 

« 2P = -1 - E ( \ 1 J S 29 " 1 ^-^- (43) 



Then the equality would, by induction, be valid for all p. 
To prove this for E 2p , consider the expression 

d 1 , d 1 

2 sinh z + cosh 2z 



dz cosh z c?z cosh z 

This is equal to zero as a result of the identity cosh 2 z = 2 cosh 2 z — 1. Expanding out sinhz and 1/ coshz and 
performing the differentiation in the 1st and 3rd terms, we have, for each power of z, 

E2 p + - +V 2 L^2Ez3l =0 (44) 

(2p-l)! + (2p-l)! + ^(2g)!(2p-2 g -l)! ' 1 > 

which can be rearranged into the form of the recursion relation (43) above. 
On the other hand we have 

P P-l P-l P-l P-l / r, 1 \ 

2p — I \ 



« 2P = -i + I2(-i) N EEE- E 

7V=2 9 2 = 1 fci = l fc 2 jv-3 = l 



2 9l - 1, 2q 2 - 1, 2fci - 1, . . . , 2fc 27v „ 3 - 1 



(45) 
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where we've separated off the N = 1 term and made the substitutions qa — >■ ki,.. . , (?2jv-i ^2^-3- We have then, 
with N = L + 1, 

a 2p = -1 + 2^ 2^ 



=1 ( 2 «! - ^'(Z© - 1)! (2p - 2gi - 2g 2 + 1)! 

1-92+1 p-gi— 92+1 p—qi— 92+1 / \ 

E E - E L 2 V' 1 , ,2 + \ • <«) 



v 2*1-1,- -,2*21,-1-1 



Setting g = ^ + g 2 — 1 , this is 



a - 1 V V f 2 {2P ~ 1)1 

U hW- 1 ) (2«)!(2p-2g-l)! 

P-9 P-9 P-9 

X 



The sum over (ft equals 2 2?_1 , evaluated by expanding (1 + l) 2q — (1 — l) 2q binomially. The remaining sum over L 
and k±, . . . , fei-i is a,2 P -2q, and we again get relation (43). QED 
The sum-over-partitions form follows in a straightforward fashion: 

£ 2 p = (-ir 1 (2p-l)! Yl \ u K u ) S 2 P -1^2m-l) km 



0<k 1 ,...,k p <2p-l 



k\ , . . . , kn 



The sign of each term in this sum is (— l)J» -1 +E) m *»'» ) while the sign in the sum-over-compositions form is (— 1) N . To 
check that the signs agree, note that 2p — 1 = J^(2m — l)fe m = 2^mk m — K. But if is the number of factorials 
in the denominator of each term, which in the sum-over-compositions form is 27V — 1. So ^ mk m = N + p — 1 and 

(_l)p-i+E mfe m = (-1)^. 

As an example of the even and the odd expansions for E n , we have: 



E w = 10! 



1 2 2 3 3 4 1 



10! 2!8! 4!6! 2! 2 6! 2!4! 2 2! 3 4! 2! 5 



/1 3 6 1 5 10 7 1\_ 

~ v _ 9! + 1^7! + 1!3!5! + 3!3 ~ TF5! ~ 1P3! 2 " + 1^3! ~ w)~ - 50 > 521 ( 49 ) 

In a similar fashion, Theorem 1 can be applied to the Euler and Bernoulli polynomials, using their generating 
functions: 

9=0 y ' U-=0 '' ) 



The results are straightforward and we omit writing out the explicit expressions. 
Bell polynomials p} are defined as 

B n< k(xx, . . . , a; n _fe_|_x) = ^fe,fe H — l-fen-it^n.Em" 1 ^ 

0< /eg , . . .,k n _k <ra 



fco!---fc p ! VI!/ "' V(»-fe + l)!/ 
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Stirling numbers of the 2nd kind are equal to the values of these polynomials at x\ = X2 = • • • = = 1: 

S{n,n-p) = B n , n _p(a;i,X2,...,Xp + i)| a . i= ... =a . p+i=1 



(n — p)! 



E 



0</ci,...,fcp<p 



n — p 
K 



K 

k\ , . . . , kp 
1 \ 101 / I 



2\ V3! 



(P+1)1 



(52) 



where we've set fc = n — p and in the last line we've summed over Uq. A similar expression for Stirling numbers of the 
1st kind can be found using their relation to the rt-th order Bernoulli numbers Q, 



s(n, n — p) 



n- 1 
P 



B 



(n) 



which have the generating function 

^' R (fe) 

E 

9 =0 



{e z - l) fe ' 



(53) 



(54) 



Using Theorem 1, we get for s(n,n — p): 



s(n, n — p) 



(n - 1) 



(n-p-1) 



1)! E 



n + JC - 1 \ i 

K ) 1 

-1 \ fel /-l 



2! 



A" 

fci , . . . , kp 

k 2 



3! 



-i 



p,E! 



(p + 1) 



(55) 



Laplace's formula (3) for the Bernoulli numbers corresponds to expression (36a) and follows from equation (2a) and 
Part (II) of Theorem 1. Other matrix representations we can derive from this theorem are: 



B 



2p 



(2p) 



2 2 p - 2 



3! 

J_ 

5! 





1 

:!! 



A 2p - (2p)! 



(2p+l)l (2p-l)! (2p-3)! 



1 



... 1 





37 



(56a) 



4! 21 



1 



(2p)! (2p-2)! (2p-4)! 
/ 1 



2! 



(56b) 



S(n,n — p) = 



(n — p)! 



n-p 



5! 1 



l l 







s(n, n — p) 



(n - 1)! 
(n — p — 1)! 



V (P+1)! pl 

/ l 

i i o 



\ 1 1 

V Ip+ttt p! 







1 \ 





(56c) 



/ 







v 



i \ 



- 1 / 



(56d) 



Note that the size of the determinant for B^ v in (56a) is (p+1) x (p + 1), (as is the one for E^ v in (56b)), compared 
with (2p + 1) x (2p + 1) using Laplace's formula. 
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IV. CONCLUSION 

We have derived compact, closed-form expressions for partition functions and for Bernoulli, Euler and Stirling 
numbers that contain only "simple" numbers and that require either finite summations or finding the determinants 
or inverses of matrices. In particular, the partition function p(n) is given directly by a sum of integers, the number of 
terms in the sum being the value at n of the pentagonal partition function. 
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Appendix A: p(5k + 4), p(7k + 5), and p(25k + 24) determinants 

For n = 5k + 4, 7k + 5, or 25k + 24, the dimension of the matrix in eq.(12) can be reduced to k + 1 by using one of 
the identities below by Ramanujan [21. Iioj: 



J2p(5k + 4)q k 

k=0 
oo 

oo 

^p(25fc + 24)<7 fe 



(9 5 ) 



5\5 



= 7 



(<Z)£, 
(g 7 ) 3 oc 



7\7 



5^-63 



k=0 



ML ' 

5 5 • 52q 



(Ala) 
(Alb) 



+ 5 10 • 6q 3 1 c12 " 4 



2 - -5--g 



- 5 7 • 63<? 2 



1!) 



where (q)oo = Y[m=i(^- ~ l" 1 )- The first identity, combined with eq.(26), gives us 



(Ale) 



p(5k + 4) 



^(-l) m (2m+ l)jm(m+l)/2 
v m>0 



('* 




E 






m= 


±1,±2,. 




1 








-6 


i 






9 


-6 


1 




10 


9 


-6 


1 


-30 


10 


9 


-6 1 


- 


-30 


10 


9 -6 1 


11 





-30 


10 9 -6 


42 


11 





-30 10 9 - 







1 \ 





■1/ 



(fc+l)x(fc+l) 



1 






-5 





(A2a) 



(fc+i)x(fc+i) 
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while the second one gives 



p(7fc + 5) 



E 

m=±l,±2,. 



^ jTn{STn— 1) /2 



+ | [ ^(-l) m (2m+l)J 7m ( m+1 >/ 2 I (/+ J2 (- l ) mj 

\ \m>0 J V m=±l,±2,... 



+ 7,j(i+ J2 (-lrj 7 "^ 3 " 1 - 1 )/ 2 ^ j> x Q 

V m=±l,±2,... / 



m(3m-l)/2 



1 \ 





■1/ 



(fe+l)x(fe+l) 



1 












1 


-8 


1 










3 


20 


-8 


1 








2 





20 


-8 


1 






8 


-70 





20 


-8 


1 




-5 


64 


-70 





20 


-8 


1 


-4 


56 


64 


-70 





20 


-8 1 


-10 





56 


64 


-70 





20 -8 


5 



(A2b) 



' (fe+l)x(fe+l) 

) 3 = l-3z + 5z 3 - 



(whcre we've used the Jacobi identity (1 — z — z 2 + z 5 + z 7 — ■ ■ ■ ) 3 = 1 — 3z + 5z 3 — 7z 6 + 9z 10 + • • • ). The matrices above 
thus consist of an LTT "base" part and a "tower" part. The coefficients of the base matrix for p(5k + 4) correspond 

to sequence A000729, the coefficients in the expansion of (IIfc(l — xk )Y ■> while the coefficients of powers of J 5 in the 
tower part is sequence A000728. Likewise, sequence A000731 gives the base matrix for p(7k + 5), while the tower part 
involves a combination of sequences A000727, A000730, and A010816. 

From the 3rd identity, the coefficients for the base matrix for p(25k + 24) are given by the expansion 

(<z) 3 ^ 



1 - 31q + 434/ - 3565/ + 18445/ - 57505/ + 70091/ + 227447/ 
(sequence A010836); the tower part is from the expansion 



3 5 

,5>\30 
oo 
,3 



63 faOgX + 5 3 • 52g + 5 s • 63/ (/^(/) c 18 



+ 5».6/(/ b oo (/)- + 5 i0 ./ (q 5 )l 

= 63 + 4988(7 + 95751g 2 + 766014/ + 3323665/ + 8359848/ 
+ 10896075/ - 6659766/ + • • • , 

and so we have, to this order, 



(A3) 



(A4) 



p(25fc + 24) 
25 



1 

-31 
434 
-3565 
18445 
-57505 
70091 
227447 



1 

-31 
434 

-3565 
18445 

-57505 



1 

-31 
434 
-3565 
18445 



1 

-31 
434 
-3565 



1 

-31 1 
434 -31 



1 



70091 -57505 18445 -3565 434 -31 



63 
4988 
95751 
766014 
3323665 
8359848 
10896075 
-6659766 



(A5) 



(fe+l)x(fe+l) 
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We list below some sample calculations of partition functions using these determinants: 



p(24) = 5- 



p(40) = 7- 



1 








1 






-6 


1 













9 


-6 


1 







= 1575; 


10 


9 


-6 


1 









-30 


10 


9 


-6 









i 










1 




-8 


1 








3 




20 


-8 


1 






2 







20 


-8 


1 




8 




-70 





20 


-8 


1 -5 




G4 


-70 





20 




8 -4 





= 37338: 



p(199) = 25- 



1 












63 


-31 


1 










4988 


434 


-31 


1 








95751 


-3565 


434 


-31 


1 






766014 


18445 


-3565 


434 


-31 


1 




3323665 


-57505 


18445 


-3565 


434 


-31 


1 


8359848 


70091 


-57505 


18445 


-3565 


434 


-31 1 


10896075 


227447 


70091 


-57505 


18445 


-3565 


434 -31 


-6659766 



3646072432125. 



Appendix B: p(5k + a) and p(25k + a) determinants 

One can ask if it's possible to fill in some of the gaps in equations (A2a), (A2b) and (A5) and to get expressions 
for p(5k + a), etc., for other values of a. In the following we will consider the problem of generalizing eqs. (A2a) and 
(A5). Ramanujan [l(| derived the relation 



(g 5 )oo (Jf + 3gJ 2 ) + gV%7f + 2qJ$) + g 2 / 5 (2J 1 2 + qjj) + q 3 / 5 ^ + g J 2 4 ) + 5g 4 / 5 
(qV 5 )oc J!-Uq + q 2 Jl 

(his eq.(20.5)), where he defined the functions Ji(q) and J2(q) by the equation 

( 9 V5) 



(9 5 )c 



= J 1 -q 1 ' s + <?' B J2. 



(Bl) 



(B2) 



Ji and J2 are series expansions in q with integer coefficients and exponents. Ramanujan then proved the identities 



Jf-llg + rt = 



(<Z 5 ) 



5^6 3 

00 



J1J2 = -1. 



(B3) 



(Ji(q) and hiq) are given by sequences A003823 and A007325, respectively.) From the first identity and from (Bl) 
we have 



5>( 



n)q 



1/5 



{Jf + 3qJ 2 ) + q 1/5 {Ji + Iqjl) + q 2/5 {2J* + q,T, 



+g 3/5 (3Ji+gJ 2 4 ) + 5g 4 / 5 



(9) 



(B4) 



and so 

oo 



(n 5 \ 5 
(Jf + 3qJ 2 ) [q )o ° 



oo 

Y,P(M+l)q k = (J! + 2qJl) 

oo 

Y J P{^ + 2)q k = (2J> + qJl) 

oo 

^p(5fc + 3) g fe = {ZJi+q4) 

fc=0 

in addition to (Ala). Let 



(<7) 6 oo 
(g 5 ) 5 oo 



(q 1/5 )oo = Gi - q 1/5 (q 5 U + ? 2/5 G 2 ; Ji • 
Then, for a = 0,1,2,3,4, 



(9 5 )c 



OO 

]>>(5fc + a) g fe = — [ F a+1 (q 5 ) a + 1 Gt a + Fi-aqiq 5 )^^ 1 ] 



where F n is the n-th Fibonacci number, with Fq = 0. We then have, making the replacement q — > J, 



p(5k + a) 



( 



(J)l + [ F a+1 {,f)^G{- a (J) + F^ a J{J^ a G a 2 +1 (J) ] x 







1 \ 




-1/ 



(fe+l)x(fe+l) 



1 












X 


-6 


1 










x 1 


9 


-6 


1 








x 2 


10 


9 


-6 


1 






x 3 


-30 


10 


9 


-6 


1 




x 4 





-30 


10 


9 


-6 


1 


x 5 


11 





-30 


10 


9 


-6 1 


x 6 


42 


11 





-30 


10 


9 -6 


x r 



(fe+l)x(fe+l) 

where the elements X n (= X^) in the tower matrix are determined for each value of a by the expansion 
Fa+^J^Gt^J) + F 4 _ a J(J 5 )^ a G a 2 +1 (J) = X I + X l J + X 2 J 2 + ■■■ 
The G's can be expressed in terms of the Ramanujan theta function: 

q 2 ,-q 3 ) 2 -"^ 2 
f(-q,-q 2 ) 

However, it is more convenient to write them as the series expansions 



f ( 2 3\2 £ ( 4\2 00 



l)/2 



Gi(g) 



-1 + V" gk(30k-l) j-j + q 2k + ^12*0+1 _ ^20fe+3 _ q 30k+7 _ (? 32fe+8 _ q 42k+U + (? 50fe+20j 



fc=0 



= 1 + q - q 3 - q 7 - q* - q 14 + q 2Q + q 29 + q 31 + q A2 - q 52 + ■ 



G 2 {q) = V" gfe(30fe+7) + _ ^lOfe+2 + ^16/0+4 + ^Ofe+ll _ ? 36fe+15 + ? 40fe+18 _ q 4Gk+23j 



k=0 



-l + q-q 2 + q i + q 11 -q 15 + q 1& -q 23 



q 37 + q i4 



q 49 + q 57 



16 



which follow directly from (B5). (G±(q) corresponds to sequence A113681, and — G2(q) to sequence A116915.) In 
column-vector form, the first terms in the expansions from eq. (B8) are: 



X (°) 



1 \ 




( 1 ^ 




( 2 ) 




( 3 \ 


1 




5 




3 




4 


9 




— 1 




5 




4 


-3 




4 




-10 




7 


-11 


; XW = 


-10 


; X< 2 > = 


3 




-16 


-10 




-7 




-9 




3 


10 




-5 




-11 




-17 


-10 




2 




-8 




-13 



V : 7 



V : 7 



V : 7 



with the coefficients for X' 4 ) being given in (A2a). 

We can generalize ec[.(A5) to an expression for p(25k + a) for a = 4, 9, 14 and 19 following Ramanujan's derivation 
of the identity (Alc)[ll|. We make the replacement q — > q 1 ^ 5 in eq. (Ala) and get 



^p(5fc + 4)< 



,k/5 



k=0 



(a 1/5 ) 6 oc 



= 5 



(g) 5 oo 1 

(<Z 5 )L (Ji-9 1 / 5 + 9 2 / 5 J 2 ) 6 ' 



To simplify the notation, we define x = qj J\ = —qJ^- Then 



Ji 



,2/5 



<h = </i(l — X 



1/5 



) = Ml + x^/^l-^x 1 ^) 



1/5-1 



where = (v5+ l)/2, the golden ratio. We have that 
1 1 - ax 1 ^ + a 2 x 2 / 5 - a 3 x 3 / 5 + a 4 x 4 / 5 



1 + ax 1 / 5 
Then 

1 



1 + a 5 x 



1 



(B10) 



(Bll) 



(B12) 



; l/5 



,2/5 



(1 + X 1 / 5 /(/))(l - f^/S) 

1 _ ^1/5 /0 + x 2/5 [J? _ x 3/5 [J? + ^4/5^,4 x + + ^2/5 + ^3/5 + ^4/5 

1 + x/(f> 5 1 — (j) 5 x 

1 - 3a; + (1 + 2x)x^ 5 + (2 - x)x 2 / 5 + (3 + x)x 3 / 5 + 5a; 4 / 5 



1- 11a; 



The denominator in (B13) is 
1 — llx — x 



jf-llq + q 2 Jl 



J! 



1 (g) 6 oo 



(B13) 



(B14) 



The numerator in (B13) is to be raised to the sixth power in eq.(B10). We define the functions H n (x) to be series 
expansions in x with integer coefficients and exponents such that 



1 - 3x + (1 + 2x)x 1/5 + (2 - x)x 2/b + (3 + x)x 3/b + 5a; 



2/5 



r 3/5 



r 4/5 



= ffi(a;) + tf2(:r> 1/5 + ff3(:E)z 2/5 
+ ff 4 (:z;):r 3/5 +#5(2^ 4/5 

Expanding the left side of this equation and collecting terms, we get 



(B15) 



H 2 (x 
H 3 (x 
H A (x 
H 5 (x 



-98a; 9 + 9939a; 8 - 107712a; 7 + 167031a; 6 - 27918a; 5 + 127011a; 4 + 160552a; 3 + 32784a; 2 + 858a; + 1 
27a; 9 - 4806a; 8 + 78758a; 7 - 171984a; 6 + 98667a; 5 + 78986a; 4 + 176592a; 3 + 52644a; 2 + 2138a; + 6 
-6a; 9 + 2138a; 8 - 52644a; 7 + 176592a; 6 - 78986a; 5 + 98667a; 4 + 171984a; 3 + 78758a; 2 + 4806a; + 27 



858a; 8 + 32784a; 7 - 160552a; 6 + 127011a; 5 + 27918a; 4 + 167031a; 3 



107712a; 2 + 9939a; + 98 



315a; s - 18640a;' + 139305a; b - 127020a; + 106425a; 4 + 127020a; 3 + 139305a^ + 18640a; + 315 
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Then, changing back to the variable q, eq.(BlO) is 



^p(5fc + 4) g fe / 5 = 

fe=0 

and we have 

oo 

]Tp(25fc + 4)g fc = 
k=0 



(<7 5 ) 



5^30 

oo 



)oo 



(5); 

5^30 



(j?Hi{q) + j! 3 H 2 (q)q^ + J 22 H 3 (q)q^ 5 + J 21 H 4 (q)q 3/5 + J? H 6 (q)qV 6 ) (B16) 



(q) 



31 



[ J 24 + 858 qJl 9 + 32784 q 2 J 14 + 160552 q 3 jf + 127011 g 4 J 4 
+ 27918 q 5 J 2 + 167031 g 6 jf + 107712 g 7 J 2 n + 9939 q 8 Jf + 98 g 9 j| 



^p(25fc + 9)g fe - 5 [ 6JP + 2138 gJ 18 + 52644 g 2 J 13 + 176592 g 3 jf + 78986 g*jf 



5\30 



fc=0 



(q) 



+ 98667 g 5 J 2 2 + 171984 q 3 JJ, + 78758 q'J 2 L2 + 4806 g 8 J 2 iY + 27 g 9 Jl 



-6 7 7 



7 t12 



„8 7 17 



„9 7221 



^p(25fc + 14)g fe = 5 [ 27 Jf + 4806 qj\ 7 + 78758 q 2 J{ 2 + 171984 g 3 j\ + 98667 g 4 J? 



5\30 



fc=0 



(q) 



78986 g 5 J 2 3 + 176592 q° J| + 52644 g Y J 2 13 + 2138 g 8 J 2 18 + 6 g 9 J 2 



,6 7 8 



7 7 13 



„8 7I8 



7231 



^p(25k + 19)g fe = 5 y-^f [ 98 J 21 + 9939 gJ 16 + 107712 g 2 j\ x + 167031 g 3 + 27918 g 4 Ji 



5\30 



fc=0 



(9) 



+ 127011 g s J 2 4 + 160552 g D J| + 32784 q' J 2 14 + 858 g 8 J 2 19 + g 9 J 2 241 

[ 63J 20 + 3728 gJ 15 + 27861 g 2 J 10 + 25404 g 3 j\ + 21285 g 4 
+25404 g 5 J 2 5 + 27861 g 6 J 2 10 + 3728 g 7 J 2 15 + 63 g 8 J 2 20 ] 



]Tp(25fc + 24)g fe = 5 : 



(<z 5 ) 



5\30 



k=0 



The expression on the right in the last equation above reduces to Ramanujan's result in (Ale) upon the substitutions 
Jl + q 2 4 = X + llq, 



J 10 + g 4 J 2 10 = A 2 + 22gA + 123g 2 , 
X 

X 4 + 44 qX 3 + 730 g 2 AT 2 + 5412 q 3 X + 15127 g 4 , 



J 15 + g 6 J 2 15 



X 3 + 33 qX 2 + 366 q 2 X + 1364 g 3 , 



J 20 + g 6 J|° 
where 

(<Z) 6 oc 



X = 



(q 



5\6 



As before, the Z n coefficients in the p(25fc + a) determinant 



p(25fc + a) = 5 ■ 



1 

-31 
434 
-3565 
18445 



1 

-31 
434 
-3565 



1 

-31 
434 



1 

-31 1 

-57505 18445 -3565 434 -31 1 

70091 -57505 18445 -3565 434 -31 1 

227447 70091 -57505 18445 -3565 434 -31 



Z 
Zi 

z 2 
z 3 
z± 
z 5 
z 6 
z 7 



(B17) 



(fe+l)x(fe+l) 
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are obtained by an expansion in powers of q of the numerators on the RHS 's of these generating- function equations. 



We have, for a = 4, 9, 14, 19, 



1 

882 
49362 
768246 
5380497 
20802996 
47413915 
46923084 



Z (4) 



V 



/ 



z< 9 > - 



/ 6 \ 

2276 
92646 
1198566 
7354172 
25710039 
51224670 
39450895 



V 



/ 



Z (14) 



/ 27 \ 

5400 
166697 
1811682 
9871992 
30828786 
55015749 
20079168 



V 



/ 



Z (19) 



/ 98 \ 
11997 
287316 
2672825 
12906450 
36553962 
54917174 
2443563 



V 



/ 



Appendix C: Y2p(n)x n determinants 

Expression (11) for p(n) can be used to express finite sums of the form ^p{n)x n as determinants. We have 



J2 P(n)x n =x k J2 x n ~ k p{n) = x k 



n=0 



n=0 



dk-i d k - 2 
dk dk-i ■ 

The tower part of the matrix can be expressed as 



1 

-1 1 

-1 -1 1 

-1 -1 



(I-J/x)- 1 







1 \ 







-1 



1 

1/x 
1/x 2 

l/x 3 

1 Xjx*- 1 
-1 \jx k 



(fc+l)x(fc+l) 



(CI) 



(fe+l)x(fe+l) 
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We now multiply the determinant in eq. (CI) by the determinant of (I — J/x), (which is equal to 1), and get 



n=Q 



n)x = x 



-1 - l/x 


1 








o 


-l + l/x 


-1-1/x 1 











l/x 


-l + l/x -1-1/x 









: 


' k -i - d k _ 2 /x 


4-2 - d k _ 3 /x 




1 





d,i, — d,i. i/x 


d.l, 1 — n /x • • • 




-1 - 


l/x 


() 


X + 1 


— X 










X - 1 


X +1 —x 










-1 


x — 1 X + 1 













-1 x-1 










2 - xdk-i d k _ 


3 - XG? fe _2 


X + 1 


—x 






_i - xd k d k . 


. 2 - xd k -i 


x-1 


X + 1 


kxk 





(fe+l)x(fc+l) 



(C2) 



d k - 
d k 

where in the last line we've expanded the (fc + l)-dimensional determinant by minors along the final column and then 
taken the factor (— l) fe x fe inside the resulting fc-dimensional determinant, multiplying each of the columns by — x. The 
final result can be written in a more compact notation as 



k 

E 

n=0 



p(n)x n = det 



q m <k+l 

-x.J T + I+ J2 (-!)' 

m>0 



jm(3m-l)/2 , jm(3rn+l)/2 _ x j(rn- 1) (3m+2)/2 _ x j(rn+l) (3rn-2)/2 



(C3) 



k x k 
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